Adiabatic quantum pumping at the Josephson frequency 
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We analyze theoretically adiabatic quantum pumping through a normal conductor that couples the 
normal regions of two superconductor/normal metal/superconductor Josephson junctions. By using 
the phases of the superconducting order parameter in the superconducting contacts as pumping 
parameters, we demonstrate that a non zero pumped charge can flow through the device. The 
device exploits the evolution of the superconducting phases due to the ac Josephson effect, and 
can therefore be operated at very high frequency, resulting in a pumped current as large as a few 
nanoAmperes. The experimental relevance of our calculations is discussed. 
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In a mesoscopic conductor in which electrons move 
phase coherently, a direct current can flow in response 
to a slowly varying periodic perturbation and in the 
absence of any applied bias. This phenomenon, known 
as quantum pumping, was first noticed by Thouless[l[, 
who analyzed theoretically the response of an electron 
system to a "traveling" periodic potential U{x — vt). 
The occurrence of quantum pumping requires that 
the periodic perturbation consists of at least two 
independent oscillating parameters X\{t) and ^(t), 
and that the trajectory representing the perturbation 
in the parameter space (XijX^) encloses a finite 
area0, Q. Indeed, the proposal of Thouless satisfies 
these requirements since even the simplest travelling 
periodic potential U(x — vt) = Uq sm(x — vt) can be 
written as U{x - vt) = Xr(i) sin(^p) + X 2 (t) cos(2^), 
with Xi <2 (t) = A\ 2 cos(^* + Xi2 ). When the cyclic 
perturbation is slower than the electron dwell time in 
the conductor, adiabatic pumping occurs and the system 
remains in thermodynamic equilibrium. In this case, the 
pumped charge can be expressed as a function of the 
scattering matrix and of its derivatives with respect to 
the pumping parameters X\ and X2 0|. 
Attempts to investigate experimentally adiabatic quan- 
tum pumping have been made using electrostatically 
defined quantum dots in GaAs-based heterostructures 
(ij. In such a system, pumping is induced by oscillating 
voltages applied to the gate electrodes defining the dot. 
Although signatures of pumping signals may have been 
observed, the experiments arc hindered by rectification 
effects originating from parasitic coupling of the ac 
signal applied to the gates Q. Many other proposals of 
devices have been put foward in the literature, in which 
different physical quantities have been used as pumping 
parameters such as a time-varying magnetic field, the 
height of a tunnel barrier, etc.jy, 0, Q. Often, however, 
these proposals do not consider the difficulties involved 
in the experimental realization. 

Here we demonstrate theoretically the occurrence of 
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FIG. 1: a) Layout of the proposed quantum pump, consisting 
of two SNS Josephson junctions in a SQUID geometry with 
a common normal region (The dark grey areas represent the 
superconductor electrodes and light grey the normal conduc- 
tor), b) Diagramatic representation of the scattering matrix 
of the device. 



pumping in a system of electrons and holes in a metallic 
conductor coupled to superconductors, where the pump- 
ing parameters are the phases of the superconducting 
order parameters in two different superconducting con- 
tacts. This system can operate at very high frequency 
without the need of feeding microwave radiation, simply 
by exploiting the evolution of the superconducting 
phases due to the ac Josephson effect. As a consequence, 
measurable pumped currents as large as a few nanoAm- 
peres, can be expected, while avoiding spurious effects 
that affected previous experiments. 

Fig. QJ, shows a schematic representation of the cir- 
cuit that we propose. Two superconducting/normal 
metal/superconducting (SNS) Josephson junctions arc 
connected in parallel via a superconducting ring, and 
their N regions are additionally coupled by a normal 
metal bridge. Andreev reflection of electrons and 
holes takes place at each NS interface, resulting in 
a phase shift of the particle wavefunction which at 
the Fermi energy is given by ±xi,2, the phase of the 
superconductor order parameter at the two different 
superconducting contacts (the sign - is for reflection 
from hole to electron; the sign + for the reverse process). 
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Hence, the total scattering matrix (Stot) of the normal 
metal bridge connecting the left and right reservoirs 
depends on the quantities Xi = e lXl and X 2 = e lX2 . We 
want to see if a direct current can flow in the normal 
metal bridge when X\ and X 2 are used as pumping 
parameters. 

The appealing aspect of such a device is the way in 
which the pumping parameters can be driven at high 
frequency, and their relative phase controlled. Specifi- 
cally, the pumping parameters become time-dependent 
when a constant voltage Vd c is present across the SNS 
junctions (e.g., by biasing the junctions with a current 
higher than their critical current), since then Xi(t) and 

Xiif) oc e % K dc * owing to the ac Josephson effect [Ioj |. 
Similarly to what happens in superconducting quantum 
interference devices (SQUIDS), the phase difference <p 
between X\ and X 2 can be easily controlled by applying 
a magnetic flux <!> to the superconducting loop, so that 
Xi = Xi + <P, with (f = 2tt$/$ ($o = h/2e). 
To demonstrate the occurrence of pumping, we model 
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FIG. 2: Pumped charge per cycle Qp as a function of phase 
difference ip between the pumping parameters, for the case 
a — b = e = 1/2 and different values of kpL. 

the system in the simplest possible way. We confine 
ourselves to the case of a fully phase coherent system 
at T = K. The normal conductor is taken to consist 
of one channel supporting ballistic motion, and the 
separation between the Josephson junctions L. The N/S 
interfaces are all supposed to be perfectly transparent, 
i.e. the probability of Andreev reflection is unity. 
The pumped current is equal to the charge pumped per 
cycle, multiplied by the pumping frequency. The calcula- 
tion of the charge pumped per cycle follows the approach 
developed by Brouwer[3[, modified to take into account 
the presence of the superconducting electrodes § . The 
relation between the charge Qp t , n pumped in one of the 
two reservoirs (labeled by m = 1,2) and the scattering 
matrix reads: 



f (dn m dX\ dn m dX 2 \ 

Qp < m = e L dt {dx- 1 ^r + dx 2 -^r)' (1) 

in which: 

dn m 1 ^ d{S t ot)i 3 , c s* , s 
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where r is the period of one pumping cycle (r = 
with wj = 2e ^ dc ). In Eq. [2j the sum over i extends 
to the electron and hole (e, h) channels in both leads. 
The sum over j is performed over the electron and hole 
channels only in the lead connected to the reservoir 
m for which the pumped charge is calculated. The 
function jij is equal to +1 when the element Si j of 
the scattering matrix corresponds to a process in which 
a current is pumped from lead m, jij = —1 when a 
current is pumped into lead m. This difference in sign 
is due to the fact that electrons and holes contribute 
oppositely to the pumped charge. Since the electron 
and hole contributions to the pumped charge could 
exactly compensate each other, it is not obvious a priori 
whether a net charge can be pumped. 
The problem of computing the pumped charge is then 
reduced to the calculation of the total scattering matrix 
Stot of electrons and holes in the normal conductor 
bridge (see Fig. [T]]), as a function of the parameters 
X\ and X 2 . The calculation is lengthy but concep- 
tually straightforward (calculations were done using 
Mathematical ) . We consider a perfectly symmetric 
configuration with two identical SNS junctions, which 
are also identically coupled to the normal metal bridge. 
For each junction, the coupling is described by a 
"beam-splitter" [llj, whose scattering matrix (£1,2) is 
assumed to be energy independent (i.e., it is the same 
for electrons and holes). We have chosen the simplest 
expression compatible with unitarity and time reversal 
symmetry. The expression reads: 




where e varies between and 1/2 (e/2 is the probability 
for an incoming particle to be deflected towards one of 
the superconductors). The amplitudes for backscattering 
a and direct transmission b across the beam splitter 
satisfy the relations a 2 + b 2 + e = 1 and e/2ab = 1, 
imposed by unitarity. For every fixed value of e two 
solutions, with a > b and b > a, are possible and we 
considered both cases (for e = 1/2 the two solutions 
coincide and a = b = 1/2). 
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Mixing of electrons and holes only occurs at the inter- 
face with the superconductors. Having assumed perfect 
transparency at the NS interfaces, the matrix describing 
Andreev reflection in the "vertical" branches of the 
circuit (see FigfTJ) depends only on the phase x 01 the 
superconducting order parameter. It reads: 
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To calculate the total scattering matrix of the device we 
first calculate the scattering matrix associated to trans- 
port across only one SNS junction. We then consider the 
two SNS junction connected in scries, i.e. we consider 
all the multiple reflection processes in the normal metal 
bridge, taking into account the corresponding dynamical 
phases acquired by electrons and holes. The result is the 
scattering matrix Stot{X\, X2) that mixes the electron 
and hole channels in reflection and transmission. 

Having determined Stot we obtain the pumped charge 
from Eqs. Q] and [2] As shown in Fig. [2] we find that, 
unless kpL = nir (with n integer), the pumped charge is 
a non-zero, anti-symmetric, and 27r-pcriodic function of 
ip, as expected. The 2tt periodicity in conjunction with 
the antisymmetry imply that the pumped charge has to 
vanish when ip — ±7r. This is the case since for tp = n 
the trajectory in the space of the pumping parameters 
(Xi, X2) does not enclose a finite area. In addition, the 
antisymmetry of Qp with respect to tp also implies that 
the sign of the pumped current changes when reversing 
the sign of the relative phase of the two superconduct- 
ing junctions. This results in the antisymmetry of the 
pumped current versus applied magnetic flux $ (that de- 
termines the phase (p), and provides a distinctive feature 
that should facilitate the experimental identification of 
the phenomenon. 

Fig. B, and C summarize the outcome of our calcu- 
lations for the different cases a > b, a = fo, and a < b. 
We first discuss the features of the results that are com- 
mon to all three cases. We always find that the pumped 
charge does not depend on the separation W between 
the beam splitters and the superconducting leads (see 
Fig. [3]). This is due to the phase conjugation [l2| of 
electrons and holes at the Fermi energy, since the dy- 
namical phase acquired by an electron propagating from 
the beam splitter to the superconducting interface is ex- 
actly compensated by the phase acquired by the Andreev 
reflected hole. In all cases the dependence of Qp on L 
is periodic for all values of ip and e, with period given 
by kpL = 7r (kp is the Fermi wave vector). This implies 
that the pumped charge is sensitive to the geometry of 
the device on the scale of the Fermi wavelength A p , indi- 
cating that charge pumping in the device considered here 
is a sample specific phenomenon. That this should be so 
is not obvious a priori : owing to phase conjugation, one 
may have expected the pumped charge to show a com- 
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FIG. 3: Color scale plots of Qp as a function of kpL and <p 
for three distinct cases of amplitude scattering on the beam 
splitter: e = 1/8 for a > b (top), a = b = e = 1/2 (center) 
and e = 1/8 for b > a (bottom). The dashed lines in the 
center panel correspond to the curves shown in Fig. 2. 



jonent independent of the precise geometry of the device 
3- 

The magnitude of the calculated pumped charge strongly 
depends on e. The maximum pumped charge is approx- 
imately 0.1 electron per cycle, for e = 1/2 (a = b). For 
small e, the magnitude of Qp decreases with decreasing e 
(and eventually vanishes for e = 0) both when a > b and 
b > a. The dependence of Qp on <p and kpL, however, 
is different in the two cases. 

When a > b and for small values of e (i.e. a ~ 1) the 
bridge connecting the two SNS junctions is only weakly 
coupled to the reservoirs, because backscattering at the 
beam-splitters is the dominant process (see Fig. [3JA). In 
this regime, sharply defined resonances appear in the con- 
ductance of the system when kpL = mr (with n integer), 
due to the presence of quasi-bound states in the bridge 
connecting the two SNS junctions. When kpL = mr the 
energy of a quasi-bound state aligns with the Fermi levels 
in the reservoirs. Interestingly, the pumped charge is also 
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significantly different from zero only when kpL is close to 
being a multiple of it. This suggests a close link between 
pumping and the presence of resonances due to quasi- 
bound states in the system, as already noted by others 
in different contexts [7[. This link is further supported by 
observing that increasing e from to 1/2 -corresponding 
to increasing the broadening of the quasi-bound states- 
results in a broader range of values of L for which charge 
pumping is observed (Fig. [3j3). 

In the case b > a, the behavior of the pumped charge 
for small values of e is qualitatively different (see Fig. 
|3]C). In this regime, the dominant process at the beam 
splitters is direct transmission. Therefore electrons and 
holes have only a small probability to be deflected from 
the normal bridge to the N/S interfaces. However, if 
they are deflected, they perform many Andreev reflec- 
tions in one of the SNS junctions before they can escape 
again to the normal metal bridge. As a consequence, 
along the dominant trajectories responsible for pumping, 
electrons and holes have a large probability to acquire a 
phase e lN% (with different, and even large, integer values 
of N), rather than simply e %x . This causes the phase de- 
pendence of the pumping signal to be richer in harmonics 
and, consequently, to exhibit very strong deviations from 
a simple sine dependence, as seen from Fig. |3p. 
Having established the occurrence of adiabatic quantum 
pumping, we briefly discuss some of the advantages of 
the proposed device. The use of the ac Josephson effect 
to generate the time dependence of the pumping param- 
eters (the superconducting phases in our case) should al- 
low operation at frequencies of several hundreds GHz. In 
fact, with superconductors such as Nb, NbN, or NbTiN, 
values for the superconducting gap A corresponding to 
frequencies in excess of 1 THz are possible, so that our su- 
perconducting pump can operate at a few hundreds GHz 
when the voltages applied across the SNS junctions is 
still sufficiently lower than A. At a Josephson frequency 
of 100 GHz, the pumped current can exceed 1 nA, which 
is easily measurable. Note that, since the pumping pa- 
rameters are coupled to the electron-hole wave functions 
via Andreev reflection, the coupling will remain good at 
these high frequencies. In addition, the fact that no ex- 
ternal microwave signals need to be fed into the circuit to 
drive the pumping parameters implies that only a neg- 
ligible high-frequency power will be irradiated, thereby 
minimizing the possibility of rectification effects known 
to cause problems in other systems [f|. For the practi- 
cal realization of the proposed superconducting pump we 
suggest the use of a ballistic InAs-based two-dimensional 
electron gas as normal conductor. Present technology 
enables the reduction of the number of conducting chan- 
nels to ps 10 [TH, which is important since the predicted 
effect is of the order of one channel. The use of InAs also 
enables the realization of the needed highly transparent 
contacts to superconductors (IBj . Furthermore, in ballis- 
tic devices in which the distance between the two SNS 



junction is L ~ 1/xm, the typical propagation time in the 
device will be of the order of L/vp 10~ 12 s (vf — 10 6 
m/s is typically realized in InAs heterostructures). This 
is ten times faster than the period of an ac pumping sig- 
nal oscillating at 100 GHz, ensuring that the dwell time 
of electrons is much shorter than the period of the ac 
pumping signal, as it is needed for the device to operate 
in the adiabatic regime. 
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